Entanglement genesis under continuous parity measurement 
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We examine the stochastic dynamics of entanglement for an uncoupled two-qubit system, under- 
going continuous parity measurement. Starting with a fully mixed state, the entanglement is zero 
for a finite amount of time, when it is suddenly created, which we refer to as entanglement genesis. 
There can be further entanglement sudden death/birth events culminating in the formation of a 
fully entangled state. We present numerical investigations of this effect together with statistics of 
the entanglement genesis time in the weak measurement limit as well as the quantum Zeno limit. 
An analytic treatment of the physics is presented, aided by the derivation of a simple concurrence 
equation for Bell basis X-states. The probability of entanglement border crossing and mean first 
passage times are calculated for the case of measurement dynamics alone. We find that states with 
almost the same probability of entanglement border crossing can have very different average crossing 
times. Our results provide insights on the optimization of entanglement generation by measurement. 

PACS numbers: 03.65.Ta,73.23.-b,03.65.Yz 



Entanglement is arguably the most fascinating aspect 
of quantum mechanics and plays a central role in quan- 
tum information science. The peculiar non-classical cor- 
relations which entanglement characterizes are the basis 
for exciting applications such as telcportation, quantum 
encryption, and many others [l| . Since the concept of en- 
tanglement appeared, there has been much research into 
its properties. One aspect of this research has been how 
to quantify the amount of entanglement contained in a 
quantum state. Many different entanglement measures 
have been introduced with varying degrees of success ■ 
In the two qubit case, the concurrence measure of entan- 
glement can be explicitly given in general for any mixed 
state 0] , making this system ideal for further research. 

In parallel with these measures, investigations into the 
dynamics of entanglement began and continue still. If 
entanglement is to be used as a resource, then its dy- 
namical evolution must be understood. Of particular in- 
terest is how entanglement decays when the entangled 
systems are coupled to the environment. It was found 
that entanglement typically decreases at an exponential 
rate, faster than single-qubit decoherence [1, H, [y, 01 . An 
important step was the realization by Jakobczyk [8| and 
Yu and Eberly d that entanglement can reach zero in 
a finite time. This behavior was dubbed "entanglement 
sudden death" [§, [I^, [HI, E [H, [H, [H and has received 
much recent attention. The ability of entanglement to 
have this non-analytic behavior can be traced back to its 
definition: two particles are considered entangled if their 
combined density matrix cannot be decomposed into a 
weighted sum of pure, separable density matrices. If such 
a decomposition can be found, then it is defined to be 
seperable, having zero entanglement. While it is impor- 
tant to understand the disappearance of entanglement, 
we must also understand how it is generated in the first 
place. This is the topic discussed here. 

Entanglement is often created via interaction-based 
proposals or with two-particle informational ambiguity 
such as parametric photon down-conversion [l6| (see e.g. 
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FIG. 1: a) Representation of the different Parity states, b) 
Diagram of a Parity Meter with two parallel double quantum 
dots on either side of a Quantum Point Contact (QPC) 



Refs. [l7|, [iSl, |l9|, |20[ for solid-state examples). From a 
fundamental perspective a fascinating alternative is the 
case in which entanglement is generated through a mea- 
surement itself either in a solid-state 
quantum optics based environment [21 
simplest case we can consider a detector which measures 
in the two-qubit Bell basis. If the system begins in a non- 
entangled state and this detector measures the system, 
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projecting into one of its basis states, then the system 
becomes entangled due to the measurement. However, 
rarely in actual experiments are the measurements per- 
fectly projective, which leads us to the study of continu- 
ous measurement. 

The subject of the creation of entanglement has been 
examined in recent years. Some papers investigate the 
idea of entanglement revivals 0, [3, EB] where the sys- 
tem is initially entangled, loses its entanglement, and 
then regains it through a revival or "entanglement sudden 
birth" event. While these effects are dynamically inter- 
esting, they might be better characterized as "entangle- 
ment sudden transfer" . This is because entanglement ex- 
ists in other degrees of freedom that becomes transfered 
to the part of the system being examined, as pointed 
out in Ref. [l^. In contrast, the current paper looks at 
measurement-induced genesis of entanglement: there is 
no reservoir of entanglement, it is creation ex nihilo by 
the act of observation. This work is closely related to 
quantum optics investigations of entanglement creation 
through a monitored decay process, see Refs. [1, H^. 

In this paper, we focus on the parity meter (a detec- 
tor which can only differentiate between qubits being 
aligned or anti-aligned) to accomplish entanglement gen- 
eration. It has been shown that parity mete rs p rovide a 
new realization of quantum computing [sO, [111 and are 
capable of building up many different types of entangled 
states [12, d^, [12] ■ Parity meters have also been recently 
been used to explore rapid two-party purification and en- 
hanced entanglement generation using feedback control 

[H. 

We explore entanglement dynamics arising as a result 
of continuous parity measurements and more specifically 
analyze the entanglement genesis event. Although our re- 
sults apply to any kind of two-qubit system, for concrete- 
ness we focus on a solid state system: a pair of quantum 
dot charge qubits coupled to a parity charge-sensor. The 
parity in this case is defined as being even (odd) if the 
qubit electrons are aligned (anti-aligned) as in Fig. la. 
There are several explicit examples of a charge-sensing 
parity meter, such as a quantum point contact with par- 
allel double quantum dots on either side (see Fig. lb) 
[13. See Refs. [H, [U, [M [11 for other physical imple- 
mentations in both charge and spin quantum dot qubits. 
The partial collapse of the system arising from the weak 
continuous measurement has recently been tested exper- 
imentally in a solid state environment [36l . [37j showing 
that the techniques required for such a measurement are 
attainable with current technology and methods. 

The results given in the text provide predictions about 
the stochastic production of entanglement from a con- 
tinuous parity measurement. This kind of entangle- 
ment generation has the counterintuitive advantage that 
stronger coupling both preserves and stabilizes the en- 
tangled state. Like the Zeno effect, any tendency for the 
entanglement to decay from coupling to other degrees of 
freedom is suppressed by the act of continuous parity ob- 
servation. Our results are also useful for the purposes of 



generating and optimizing entanglement production. For 
example, it may be easy to prepare certain initial states 
that minimize the time to cross the entanglement border 
(|27p . given a fixed probability of entanglement genesis 
(I24I251) . 



I. MODEL 

Under the assumption that the qubits are spatially 
separated and screened by the quantum point con- 
tact (QPC), any direct qubit-qubit interactions are 
suppressed leading to the qubit Hamiltonian H ^ 
[eial + 62(7^ + Aio-^ + A2(tI] /2, where Aj is the tun- 
nel coupling energy between the dots, which is set to be 
symmetric (Ai = A2 = A). For simplicity assume that 
the qubits are tuned such that the energy asymmetry 
between the levels (e^) is set to zero. This Hamiltonian 
gives a qubit time scale, Tq = 27r/A, and has four energy 
eigenvectors (-01, i/'2, "03, i/'4), 
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(|ll) + |10) + |01) + |00))/2, 
(|ll)-|10)-|01) + |00))/2. 



where |10) refers to the first qubit being in the "up" state 
and the second being in the "down" state (etc.). The first 
two solutions share the interesting property that they 
both have a definite parity and therefore are stable under 
the measurement dynamics as well as the qubit Hamil- 
tonians. The second two are stationary under Hamilto- 
nian dynamics, but not under measurement dynamics, 
since they are not parity eigenstates. Another interest- 
ing aspect about the first two states is that these are 
Bell states. Therefore, by combining the dynamics from 
these two processes (measurement and Hamiltonian evo- 
lution) we end up with two stationary states that are 
both maximally entangled and a third possibility, where 
'03, '04 continually rotate among themselves and are sub- 
jected to a continuous combination of noisy measurement 
and Hamiltonian dynamics. A time-averaged electrical 
current alone can differentiate between these three be- 
haviors of the system at sufficiently long times [2^ . 

The coupling between the qubits and the QPC can be 
modeled with the interaction Hamiltonian, 



i/,„t = {^E/2)alal, 



(2) 



where Ai? is a charging energy operator, whose expecta- 
tion value is equal to the difference in charging energy of 
the even and odd parity classes. In the tunneling regime 
the operator AE can be associated with the standard 
input variables Aaj^a^ + H.c, where A is the coupling 
strength and ^ and ajij^ are the creation and annihi- 
lation operators associated with the right or left side of 
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the tunnel junction. For more details on this model of 
the interaction see (2^ . 

For the systems considered here the output signal of 
the parity meter is the electrical current, I{t), through 
the quantum point contact. This current is modeled fol- 
lowing Korotkov [lil by I{t) = J2kPkkIk + ^(t), with 
Pkk being the density matrix clement associated with the 
state IV'fc) and ^{t) being the random white shot noise of 
the QPC, described with a normal distribution with spec- 
tral density, 5*0 defined as {^{t)^{0)) = 2SoS{t). Each I, 
corresponds to the average current that the correspond- 
ing state would produce. Defining the current pro- 
duced from an even (odd) parity state as Ie.o we can 
introduce / = {Ie + Io)/'^ and A/ ^ Ie — lo- This al- 
lows the definition of dimensionless shifted currents with 
lit) — / = x{t)^, where x = ±1 coincides with the 
current from an even/odd parity state. From this model 
an important time scale arises, the measurement time: 
Tm = 45*0 /A/^. This is the typical time needed to differ- 
entiate between the noise and the signal from the system. 
With the measurement time and the qubit time a useful 
parameter can be defined as K = Tq/TM, which can be 
varied by tuning parameters of the system. K can be 
increased through a stronger coupling between the sys- 
tem and detector, as well as a higher potential barrier 
between the quantum dots themselves. 

To analyze the entanglement we follow Wootters Q 
and define concurrence as C = max[0,A], where A = 
•s/AT — -I- a/As -I- v^Al) and corresponds to the 

eigenvalues of the matrix p{ay ® ay)p*{ay ®(Ty), ordered 
such that Ai > A2 > A3 > A4. Previous works have 
achieved significant analytical understanding by applying 
concurrence to a density matrix that maintains a general 
X form in the computational basis (|00), |lf). |01), and 
|10», 



to formulate a simple expression for the concurrence, 
C = 2max[0, \p2z\ - ^piipa, \pii\ - ^/p^UpE] [H- Al- 
though this expression is very useful, the density matrix 
we presently consider is not trapped in an "X" state in 
the computational basis. However, using the Bell Basis 
defined as: 



K) = (|ll)-|00))/^/2, 

K) = (|ll) + |00))/V2, 

K> = (|10) + |01))/^/2, 

\u,) = (|10) - |01))/^/2, 



(4) 



provides a similar situation. The isolated qubit Hamilto- 
nian expressed in this basis is given by 



H 




(5) 



The structure of this Hamiltonian together with the fact 
that the first and last Bell states share a definite parity 
guarentees that if the system starts in a Bell basis X-state 
^ , then it will remain in one. This fact again simplifies 
the calculation of the concurrence. This is an important 
difference between previous work and our current inves- 
tigation which does provide significant differences. 
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The eigenvalues for p{(jy®ay)p*{(jy®(jy) corresponding 
(3) to this general X-state are 



'K,b = 2 V (Pll + Pii + Pli + PWiPll + Pii ~ Pl4 - Pm) T V {Pll - Pa + Pli - Pl4)(Pll - P44 - Pli + P*u) 



{P22 + P33 + P23 + P2z){P22 + P33 " P23 " P23) T V (^22 - P33 + P23 " P23)(P22 - P33 ^ P23 + P23) 

I 



(6) 



From these eigenvalues it is clear that ^/Ab > V%i and 
^/Xd > •\/A^, making either -y/Ah or compete for 

the role of \/Xi. With these two orderings we have 
two possible A functions that compete to be the con- 
currence. In spite of the lengthy expressions for the 
•\/A's, a pleasing simplicity emerges to give a concise 
and physically appealing result. We expect that this 
result will be of interest in its own right. We assume 
that the initial state is fully mixed (p = 1/4), which 



means initially pi4 and P23 = and from the dynamics 
([5]) it can be seen that puit) = for all t and that p23 
will remain completely imaginary. With this assumption 
one possible A function is split further due to the term, 
VAPii^^°44)^ = |pii— /544I, yielding three possible values 
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for A: 



Ai = 2pn - 1 
A2 = 2p44 - 1 



A3 = V(/022 -P33)^ +4|p23p + P22 + P33 - 1- 

The correct A; to use for the concurrence 

C = max[0, A] 
is the maximum of the three 

A = max[Ai, A2, A3]. 



(7) 



(8) 



(9) 



Physically these three different functions can be under- 
stood as the concurrence coming from the three different 
possible asymptotic outcomes of the parity meter, either 
one of the stable solutions with even parity (Ai) or odd 
parity (A2), or the oscillatory solution (A3). In the fully 
diagonal case we have the simple result A = 2 max pa — l. 

For further examination of the measurement dynamics, 
we follow Korotkov, Ruskov, and Mizel [H,!!^ and write 
a stochastic differential equation (in Ito representation) 
governing the evolution of the system: 



So 



45o 



if"'' 



(10) 

where again I{t) = Y^kP^klk + ^{t), with ^(t) being 
the random noise from the measuring device and /; cor- 
responds to the current signal that the corresponding 
state \ui) would produce (in our normalized units the 
even states Ie ~ Ii ^ h = ^ and the odd states 
Iq = I3 = li = —I)- The rates 7^ are the environment- 
induced decoherence rates, which we neglect in order 
to focus on the fundamental physics. Due to the non- 
analytic nature of the concurrence we cannot create a 
single differential equation to govern its evolution, so we 
will initially concentrate on numerical analysis. 



II. NUMERICAL ANALYSIS 
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Starting from an initial fully mixed state, pin = 1/4, 
we can allow the system to evolve governed by ([TO)). A 
typical run is displayed in Fig. 2. We can see several in- 
teresting features in these runs. First, the system starts 
out unentangled due to the initial fully mixed state and 
then after some amount of time experiences an entangle- 
ment genesis event. Understanding this time scale is a 
primary goal of this investigation. Secondly, we see that 
entanglement sudden death behavior, first examined in 
optical systems Q, is a common occurrence in this 
system. This sudden death behavior can be understood 
as the system initially being projected into one of the 
three outcomes from the parity meter, and then due to 



FIG. 2: The top figure shows the evolution of the concur- 
rence throughout a typical run with K = 0.3, plotted as a 
function of time in units of the qubit time Tq. The dashed 
lines represent the bounds on concurrence (C = being un- 
entangled, C = 1 being fully entangled). There is always a 
period of quiet time punctuated by the entanglement genesis 
event. Notice that there are several occurrences of entangle- 
ment sudden death/birth. The middle figure shows the inte- 
grated detector output, with the fines representing the output 
from a definite purity state with no noise. The bottom fig- 
ure shows an average electrical current, with the dashed lines 
representing the average current from states of definite parity. 
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FIG. 3; Average A function with K = 0.03, 0.3, and 30, aver- 
aged over 1000 runs. It is useful to consider K — 0.03, 0.3, 30 
to characterize when the case of a weakly coupled system 
{K — 0.03), a moderate coupUng {K = 0.3) and a system 
that is in the strong couphng regime {K = 30). The entan- 
glement genesis time is when A crosses 0. 



FIG. 4: Histogram showing the frequency of particular en- 
tanglement genesis times from 10,000 runs, with K = 0.3. 
Each bar represents 0.2 Tq. Notice the gap in entanglement 
times at the beginning as well as the long tail where the en- 
tanglement genesis time can be larger than lOT, due to the 
competing A functions. 



noise the system will begin to drift back toward a differ- 
ent outcome causing the system to lose its entanglement 
for a period of time. In addition to the sudden death, 
we also see the recovery of entanglement through sudden 
birth transitions numerous times during the time evo- 
lution. The entanglement is always regained eventually 
since each asymptotic solution to the dynamics is a max- 
imally entangled state. 

Another quantity that we can examine with these nu- 
merical simulations is the average A function, seen in 
Fig. 3. We can again see a quiet period before entan- 
glement is generated (when A < 0). This entanglement 
genesis time is dependent on the coupling strength of 
the meter. The average concurrence asymptoticly ap- 
proaches 1 in an exponential fashion. Of interest is the 
spread around this average genesis time. Running 10,000 
realizations of the quantum measurement and plotting a 
frequency histogram (see Fig. 4) we find the following: 
The distribution is zero up until some finite point, in- 
dicating that there is always a quiet period before the 
entanglement is generated. This can be understood as 
the system having to traverse the Hilbert space from its 
originally fully mixed state and cross the entanglement 
border to the closest state which is entangled. Using the 
trace distance [l|, defined as d = iTr|/9 — cr| where p and 
cr are two states, the closest entangled state that has a 
non-zero trace distance of c? = 1/16 is 



a/4 

1/4 i/A 

-i/A 1/4 

l/4y 



(11) 



This state is on the boundary, A = 0, of entanglement 
and it should be noted that this state can be achieved 
with just the qubit Hamiltonians. However, without the 
measurement dynamics, this state will never cross the 



border and generate finite entanglement. 

One further numerical result is obtained in the in- 
teresting regime where the coupling strength between 
the measurement device and the system is very strong 
{Tq/TM ^ 1), rapidly projecting the state. We expect 
to see quantum jump behavior similar to the Quantum 
Zcno Effect (see Ref. H^). Studying the numerical re- 
sults in Fig. 5, this type of sudden shift in the system 
can be seen as well as other interesting behavior. Very 
soon after the parity meter is turned on, A leaps from 
its initial value —1/2 to zero, corresponding to the fully 
mixed state being projected quickly into the odd parity 
subspace (/ = —1). Since the system started out in a 
fully mixed state, this parity subspace has both states 
in the subspace almost equally occupied, causing the A 
function to be around A = 0. Secondly, after this initial 
leap to zero the concurrence has a period of bouncing 
around near zero. This can be understood as the qubit 
Hamiltonian rotating a small probability amplitude be- 
tween \u2) and I Ms), and then the measurement collapses 
the system back into the odd subspace. A final quantum 
leap is seen when the small rotation to the even parity 
subspace causes the system to collapse toward this new 
outcome, resulting in the jump in concurrence and leav- 
ing the system in the oscillatory solution with C = 1. 



III. ANALYTICAL APPROACH 

In order to have a better understanding of this physics 
we can take a simplified analytic approach. Let us as- 
sume that we have projective measurements giving either 
"even" or "odd" at each step, with unitary Hamiltonian 
evolution occurring between them. This will allow us to 
examine the average concurrence at each step and derive 
an expression for the concurrence after an arbitrary num- 
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FIG. 5: These figures demonstrate a run with if = 30 where 
we are in the Quantum Zeno Regime. The top figure shows 
the evolution of concurrence. The middle plot again shows the 
integrated output and the bottom plot shows the averaged 
output. The concurrence rapidly leaps to C = 1 under a 
projection into the even sub-space. 



ber of measurement steps. Starting with a fully mixed 
state p = 1/4, we have an initial concurrence, C = 0. 
The probability of measuring an even parity result is 
Pu + P22 and the probability of an odd parity result 
is P33 + P44- Since we are now dealing with projective 
measurements, each measurement will fully project the 
system into its respective parity subspace. The first mea- 




10 

Time (Tq) 

FIG. 6: The ensemble average of A is plotted versus the qubit 
time Tq. The solid line represents the numerical results from 
10,000 runs with K = 30. The dotted curve is the simplified 
projective analytic result Eq. (|16|l with S = n/K. Both curves 
start out at A = —0.5, however the projective model jumps 
to zero after one measurement time {Tq/30). 



surement will result in 



Pi5 = (1/2) diag(l, 1,0,0), 



(12) 



with probability 1/2 corresponding to the first result re- 
turning "even" or 



po = (1/2) diag(0, 0,1,1), 



(13) 



with probability 1/2 corresponding to the first result re- 
turning "odd" . Each of these results has A = 0, so 
the average concurrence is zero for the first measure- 
ment (n ~ 1). Now the qubit Hamiltonian dynamics 
comes into play and mixes the parity. Assuming a small 
time step 5t such that S = St A/2 ^ 1 and that the first 
measurement returned an even result, we can write the 
evolved density matrix as 



Pe = 7^ 










cos^ 5 
-i cos 5 sin 5 




Q\ 

i cos S sin S 






sin^ S 



(14) 



0/ 



From here on wc can just consider this case because 
if the first result were odd, the analysis is a mirror im- 
age of that below, and gives the same average. A sec- 
ond projective measurement will again collapse the sys- 
tem into a parity subspace, where there is a probability 
Po = (l/2)sin^(5 that the result be odd, preparing the 
Bell state lua) with concurrence Co = 1- However, for 
(5 <C 1 this probability is small, and most of the time 
the system will be projected back into the even subspace 
with probability pE = {1 + cos^ S)/2, preparing the state 



1 



1 + cos^ 6 



(\ 0\ 

cos^J 



Vo 0/ 



(15) 
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which has a concurrence of Ce = 2 maxj p'^ ■■ — 1 = 
(1 — cos^(5)/(l + cos^(5). This situation is a two-qubit 
analog to the Zeno effect, and similar to the discussion 
in Ref. [131 ■ Once there is a result that differs from the 
earlier outcomes of the parity measurement, the system 
is then always trapped in an oscillatory combination of 
1^2,3), from then on contributing C = 1. Therefore the 
probability of having a sequence of the same results is the 
important quantity to follow in the evolution, as well as 
the value of the concurrence associated with those states 
in the sequence. The average concurrence (C) after the 
second measurement {n = 2) is given by (C) ~ PeCe + 
poCo = 1 — cos^ 5. Generalizing this approach, we find 
that the average concurrence as a function of the number 
of measurement steps n is 



{C){n) 



1 



s2("-i) 5. 



(16) 



This solution can be compared to the continuous Zeno 
model by considering the numerical average of A over 
many realizations with K = 30 by choosing 5 — t: / K 
and comparing it with the result plotted in Fig. 6, with 
reasonable agreement. We note that the projective model 
rises to a non-negative concurrence after the first mea- 
surement, while the continuous measurement model takes 
longer to create entanglement initially, and consequently 
lags behind the projective model. 

We next consider the case when the measurements are 
weak. Unfortunately, following the same calculations 
as above for weak continuous measurements are rather 
cumbersome. To make further analytical progress, we 
will adopt a strategy similar to Korotkov and one of the 
authors in Ref. 40| to investigate entanglement border 
crossing. We proceed by ignoring the unitary evolution 
from the qubit Hamiltonian dynamics and concentrate 
on the measurement dynamics alone. We consider an 
arbitrary diagonal state in the Bell basis 



(17) 



A weak continuous measurement of the current yields 
new information about the system. As a result of this 
information, the density matrix of the system should be 
updated. Since we are neglecting the qubit Hamiltonian 
we can use the classical Bayes rule to update the density 
matrix. This is done by interpreting the combination 
Pe = Pii + P22 as the probability to be in the even sub- 
space, and the combination po = P33 + Paa as the prob- 
ability to be in the odd subspace. Given the observation 
of an electrical current /, these probabilities are updated 
using the conditional probabilities of finding the current 
/, given a definate parity, Pe,o{I): via Bayes formula, 
p'e = Pe{I)peI [Pe{I)pe + Po{I)po\, and the same for 
p'q . We take Pe o {I) to be Gaussian distributions given 

by 



(Pll 














P22 



















Vo 








Pii 



Pe. 



^t/TiSo exp -{I -lE,o?tlSo 



(18) 



where Ie,o is the average current from an even (odd) 
parity state, and introduce the rescaled measurement re- 
sult 



7(t) = (/-/o)(iA//5o), 



(19) 



to find that the updated density matrix conditioned on 
the measurement result 7 is given by 



1 

















P22e^ 














/933e-^ 














pne 



(20) 



with N = (pii + p22)c' + (P33 + P44)e Notice that 
as 7 ±cx) the system collapses into the even (odd) 



subspace as expected. From (O, A is 



A 



max[piic'',P22e'^,P33C '^,^440 ^ 



(21) 



Entanglement border crossing happens whenever this A 
expression crosses zero, which can correspond either to 
sudden death of entanglement or entanglement genesis 
depending on the direction. In order to explore these 
phenomena we need to determine what value of 7 will 
correspond to A = 0. Setting (PT|) equal to zero we can 
solve for 7, producing two solutions, 7 = ri or r2. 



rr = iln 

r2 = -iln 



P33 + P44 



2max[pii,p22] - Pll - P22 
Pll + P22 



(22) 



.2max[p33, P44] - P33 - /944. 

These two solutions correspond to the possibility that 



a) 1 



P44 

0.5 



b) 1 



P44 

0.5 



0.5 

Pas 



0.5 

'^33 



FIG. 7: Illustration of entanglement border crossing, condi- 
tioned on the measurement result 7. Initially the states in 
the lower quadrant (marked by dashed lines) are unentan- 
gled, while states outside the lower quadrant (but still in the 
physical lower triangle) are entangled, (a) 7 > 0: Some ini- 
tially entangled states experience entanglement sudden death, 
(b) 7 < 0: Some unentangled states experience entanglement 
genesis. Entangled states are shaded in purple, and the bor- 
der is given by straight lines originating from the center of 
the square. 

starting from an unentangled state, one can become en- 
tangled via pure measurement dynamics by collapsing to 
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either the even or the odd parity subspace. From Eq. ((22|) 
we can see that if pn = P22, or ^33 = P44 we will have 
an infinite ''1(2), which makes physical sense. Consider- 
ing the case pn = P22, if the system completely collapses 
into the even parity subspace (which would correspond 
to 7 +00) this results in the density matrix elements 
P'li ~ P'22 = 1/2 corresponding to an equal mixture of the 
first two Bell States with A = 0. Therefore if one starts 
with pii = P22, then the even route to entanglement is 
blocked, leaving only the odd route. 

We now analyze the measurement dynamics and find 
the probabilities and mean time needed to cross the en- 
tanglement border, A = 0. We begin by noting that the 
conditional probability distributions (|18p arc solutions to 
the two diffusion equations 



E,0, (23) 



where ve,o = ±1/7*/ correspond to the different drift 
velocities depending on the parity, and D = \/{2Tm)- 
The probability distribution of the measured current, 
P(I) = (pii + P22)Pe{I) + {P33 + P44)Po{I), docs not 
involve any off-diagonal density matrix elements, and 
therefore the calculation is identical to a classical sys- 
tem with probabilities pe,o of being in state "even" or 
"odd" . We may therefore solve the problem for each case 
separately, and then weight the solutions with the prob- 
abilities pe,Po of occupying a given parity. 

We wish to calculate the probability that the random 
variable j{t) crosses the specific values ri^2 by solving 
the diffusion equations using a Green function approach. 
This approach places absorbing boundary conditions at 
7 = ri^2 in order to find the first passage time statis- 
tics. The details of this calculation are in the Appendix. 
Several important results emerge from the analysis. As 
described previously, wc neglect singlc-qubit Hamiltonian 
dynamics (A ~ 0) and assume that the initial state is di- 
agonal in the Bell basis with pu = P22 7^ and (033 7^ (044. 
The condition pn = P22 7^ together with the normaliza- 
tion of the state reduces the number of independent state 
variables to 2, which we take to be P33, P44. States within 
the lower quadrant of the P33 — pa plane are seperablc, 
while physical states outside the quadrant are entangled. 
The entanglement border is drawn as a dashed line in 
Fig. 7. As the measurement progresses, the updated state 
(PO)) can cross the entanglement border, depending on the 
value of 7. In Fig. 7 we indicate the initial states that 
cross the entanglement border for (a) 7 > (where some 
states experience entanglement sudden death), and (b) 
7 < (where some states experience entanglement gen- 
esis). The entangled states are shaded in purple, and 
the border is given by the lines P44 = 0^33 -|- 6, where 
a = (1 ± e2^)/(l T e^T), b = 6^^/(1 =f e^T) for the cases 

P44 ^ P33- 

We now present our results for the specific case dis- 
cussed above. The probability to begin in an unentan- 
gled state (in the lower quadrant of the P33 — P44 plane) 
and to have an entanglement genesis (EG) event at any 



Unphysical 




P44 
0.5 







■ 






1 









0.5 /)33 



FIG. 8: Top: The average entanglement border crossing time 
Tc is plotted in units of the measurement time, Tm, as a 
function of P33 and P44 for states where pn = P22 7^ and 
P33 7^ P44. The contours are set at time increments of 0.2 with 
dark blue (0) to yellow (1.4) (white exceeds 1.4). Bottom: 
Probability of an entanglement border crossing {Peg or Psd 
depending on the initial state) as a function of P33 and P44. 
Notice the change in behavior between EG and SD across the 
lower quadrant boundary. The contours are set at every 0.1 
probability increment with dark blue to yellow representing 
to 1. 



point in time is given by 

Peg = 2max[p33,p44]. 



(24) 



If we start from an entangled state (outside the lower 
quadrant) we find the probability for a sudden death 
(SD) event at any point in time to be 



P' 



SD 



P, 



EG- 



~ P33 - P44 
|P33 - P44I 



(25) 
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It is remarkable that the results naturally break into two 
parts, depending on the direction of the border cross- 
ing, and further arc independent of all details except the 
initial density matrix elements. The reason for the differ- 
ence between Peg versus Psd can be intuitively under- 
stood by considering the extremes. The entangled states 
close to the line ^33 -t- P44 = 1 can only experience sud- 
den death if the system collapses toward the even parity 
state. However, P33 -I- P44 is the probability to have odd 
parity, so these states are unlikely to experience sudden 
death. The other extreme corresponds to unentangied 
states close to the origin. These states can only experi- 
ence entanglement genesis if the system collapses toward 
the odd parity state. However, being close to the origin 
means that the probability to be in the odd subspace is 
small, therefore entanglement genesis is unlikely for these 
states. 

The average time, Tcj it takes the system to cross 
the entanglement border conditioned on the fact that it 
will cross can be determined from the first passage time 
probability (calculated in the Appendix) 



EG 



fpt\c 



\r2\ 



V27rT3 



exp 



2t 



and is found to be 



In 



1 - P33 + P44 
|P33 - P44I 



(26) 



(27) 



from (jA.lip and for either sudden death or entanglement 
genesis. Plots of Tc and PeGi Psd can be seen in Fig. 8. 
The above equations relate to the system being collapsed 
into the odd subspace to create or lose entanglement, but 
by symmetry if the assumptions involved are reversed: 
P33 = P44 7^ 0, pii ^ P22, then pii^22 would exchange 
roles with ^33,44 in the above equations, and the system 
would collapse into the even parity subspace to create or 
lose entanglement. 

For clarity and further insight, we now consider the 
specific cases of: 

i) an almost fully mixed state with 

Pii = P22 = 1/4 + e, P33 = 1/4 - 3e, and P44 = 1/4 + e, 

(28) 

ii) a near-Bell state with 



Pii = P22 = f , P33 = 0, pm = 1 - 2e, 



(29) 



and a state close to the boundary on cither side (non- 
entangled) iiia) with 

Pii = P22 = 1/4, P33 = e, P44 = 1/2 - e, (30) 

and (entangled) iiib) with 

Pii = /322 = 1/4 - e, P33 = e, P44 = 1/2 -I- e. (31) 

For these states, we find that as e ^ 0: 

i) Almost fully mixed: 



1/2 + 2e 
In 



1/2, 



l/2+2e 
4e 



ii) 



iiia) 



iiib) 



Near-Bell state: 
Psd = 4e ^ 0, 



Tc 



_ Tm. 
2 



In 



l-2e 



Non-entangled state near border: 

Pbg = 1 - 2e ^ 1, 



In 



1/2 



1/2-2C 



0. 



Entangled state near border: 
Psd = (1/2 - 2e) (l + 



It is interesting that in case (i), the state will become 
entangled with 1/2 probability, but the time required 
for this to happen diverges logarithmically. The reason 
for this can be understood from Fig. 7. For the near- 
fuUy mixed state, half of the time the system collapses 
into the odd parity state, so 7 — * —00. In Fig. 7(b) 
the entanglement border is pinching in on the state lo- 
cated near (1/4,1/4). The entanglement border ap- 
proaches the state in an exponentially decaying way, and 
that is why the crossing time is logarithmically diver- 
gent. Another surprising result from these calculations 
is that some states take a long time to cross the entan- 
glement border, but have the same crossing probability 
as states that take a very short time. For example, if 
Pii = P22 = -25, P33 = .49, P44 — .01, we find that 
Tc = (rM/2)ln[50/48] w 0.02Ta/ and Pc = .98. How- 
ever, with a different state, pn = P22 ~ -02, P33 = .49, 
P44 = .47, we find the same Pc but the crossing time be- 
comes Tc = {Tm/2) ln[2] « 0.35Ta/. Two states with the 
same probability of crossing have an order of magnitude 
difference in their average crossing times. 



IV. CONCLUSION 

In this paper the entanglement genesis and sudden 
death behavior of a two-qubit system coupled to a par- 
ity meter has been modeled and investigated. We found 
from numerical simulation that both effects are generic. 
We developed a simplified analytical approach to look 
at the behavior of the concurrence as well as the typi- 
cal time scales involved in the entanglement generation. 
We found the probabilities of entanglement genesis and 
sudden death (|24l25p . as well as the average time of en- 
tanglement border crossing l|27p . When starting from a 
fully mixed state there will always be some duration of 
time before entanglement is generated due to the dis- 
tance between the initial and nearest entangled state in 
the Hilbcrt Space. This entanglement genesis time can 
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vary greatly for different states, even though the proba- 
bihty of having such an event remains comparable. The 
detailed understanding of entanglement creation in these 
systems sheds light on the fundamental nature of contin- 
uous entanglement dynamics, as well as providing impor- 
tant aids in the development of methods to control and 
optimize entanglement generation. A fascinating topic 
for further research is the use of feedback to steer the en- 
tanglement production. A primary goal of this research 
would be to prevent entanglement sudden death, or to 
shorten the entanglement genesis time. 



APPENDIX 

The purpose of this Appendix is to explain the Green 
function approach to calculating the probabilities of 
entanglement border crossing, and the associated first 
crossing time statistics. We begin by rescaling time in 
units of the measurement time, t = Tmt, so the con- 
ditional probabilities (jlSp of a measurement outcome 7 
associated with a given parity subspace are given by 



^E.O 



(7) 



1 

27rr 



exp 



(7Tt) 
2t 



(A.l) 



These distributions arc solutions to two random walks 
with drift velocities ve,o = ±1 and a diffusion constant 
of -D = 1/2. These random walks are described by the 
Fokker-Planck equations 



i = E,0. 



(A.2) 



If wc want to look at the first passage statistics relat- 
ing to A crossing zero then we must first find the Green 
functions, Gi which satisfy (jA.2p . with the initial condi- 
tion that 7 = 0. Following the argument given in the 
body of the paper, the solutions of the Green functions 
Gi{"f,T) for a definite parity i = E,0 are weighted by 
the even/odd parity probabilities pe = Pii + P22 and 
Po = P33 + Paa respectively in the final results. We 
place absorbing boundary conditions at 7 = ri and r2 
to account for the fact that we are only interested in the 
first passage time (for more details see Ref. [4l|). This 
gives us the condition that 0^(7 = ri^2iT) = 0. Fol- 
lowing the discussion in the text, we now focus on the 
case pii = P22 7^ and P33 7^ p^A in order to simplify 
the analysis, because ri —* oo giving only one finite ab- 
sorbing boundary condition at 7 = r2. We recall that 
gaussian functions are solutions to differential equations 
of this type so we can find the Green functions by exam- 
ination: 



G,(7,r) 



1 



I exp 



(7 - ViT)' 



ADt 



(A.3) 



1 — cxpl — 



r2 - r2l 
Dt 



The above equation satisfies the differential equation 
of (jA.2p as well as the necessary boundary and initial 



conditions. We now focus on the entanglement gen- 
esis case and take r2 to be negative. The function 
P<oi''') ~ Ir^ is defined as the total proba- 

bility that A is less than zero at time r for a given initial 
state. In the sudden death case, r2 is positive so the sys- 
tem is initially entangled. Here, we define the function 
P^q(t) = G'i(7,r)c?7 to be the probability for A re- 
maining larger than zero. For definitcness, wc now focus 
on the entanglement genesis (EG) case (the sudden death 
(SD) case is directly analogous). Since there is one finite 
absorbing boundary, the only place to lose probability is 
at 7 = r2. Therefore, the first-passage time probability 
distribution P|pj(t) is given by 



P}pt{r) = -drPUr), 

= -dr I G,(7,r)d7, 

) 

drGidj. 



(A.4) 



From (|A.2[) we have that 9rGj = -Vid^G^ + Dd'^Gi so 
we can rewrite the above equation as 



Pfvtir) 



v,G, + Dd^GAdj 



- v,G, + Dd^G, 

{r2 - v^tY 



DdjGi\r2 1 



(A.5) 



1^2 



: exp I 



ADt 



where we have used the boundary condition (7 = r2 ) = 
0. The SD case is accounted for with the absolute value 
of r2. Integrating this result over all positive time will 
give us the probability that the A = point is crossed 
eventually. 



exp 



r2Vi 



\r2Vi 



2D 



(A.6) 



The EG case corresponds to r2 < 0, and we can therefore 
find Peg (the probability of entanglement genesis given 
parity state i = E,0) is given by 



P, 



EG 



for i = E, 
for i — O, 



(A.7) 



while the SD case corresponds to r2 > 0, and we can 
therefore find Pg^ (the probability of entanglement gen- 
esis given parity state i = E,0) is given by 



P. 



SD 



for t 
for i 



O. 



(A., 



The probability of 1 makes physical sense because for 
the (odd, even) parity states, the (negative,positive) drift 
velocity W2 = Tl will eventually send 7 past r2. The 
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exponentially suppressed probabilities correspond to the 
noise fighting against the drift term. 

We can now weight these probabilities with the ini- 
tial density matrix elements to get the full probability 
of crossing, Pcross = PePctoss + PoPcross- Combining 
results, we find the entanglement genesis probability is 
given by 



P; 



EG 



{pii + P22)e """^ + (/O33 + P44) = 2max[p33, P44], 

(A.9) 

where we used the threshold (HH). Analogously, the sud- 
den death probability is given by 



Pi 



SD 



(Pll 
(Pll 



P22) + (/O33 + P44)e' 
2max[/933,p44] 



-2r2 



P22) 



\P33 - Pii\ 



(A.IO) 



Accounting for normalization recovers Eq. (|25p . 

We can rcnormalizc (|A.5p with (|A.6[) to get the prob- 
ability distributions Pfpt\c '-'^ crossing as a function of 
time, conditioned on the fact that they will cross. This 
will allow us to determine the first passage time statistics. 



These conditional probabilities are 



pi 



Using this probability density, the average first passage 
time is given by (t,) = dr tP'^^^^^ = |r2|/|w,|. When 
this average time is weighted by the initial density matrix 
elements, the total average time until the border crossing 
is 

(t) = Pete + POTQ = \r2\- (A. 12) 

Restoring physical time units gives the result (|27|) . 
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